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ABSTRACT 


Coherent light beams propagating through the atmosphere undergo 
considerable phase perturbations due to fluctuating temperature structures in 
the atmosphere. Understanding and measuring these structures on a real-time 
altitude dependent basis is inherent to successful deployment of ground based 
lasers and particle beams. 

One method used to detect these temperature structures is an acoustic 
profiler, or echosounder. Of immediate interest is the ability of high frequency (5 
KHz) planar array echosounders to rapidly detect low level turbulence (below 
200 meters) and quantify the results. 

This thesis involves design improvements of previously developed 
echosounder arrays and associated software. Particulary, this thesis 
demonstrates that tighter packing of elements in a planar acoustic array 
produces better side lobe reduction than less densely packed arrays. This 
results in higher energy density in the main lobe and increased performance. 

Also included in this thesis is a computer method which allows relatively 


accurate beam pattern prediction from any given planar array. 
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l. INTRODUCTION 


Optical degradation of coherent lightin the atmosphere is strongly a 
function of small scale temperature gradients which change the local index of 
refraction in the atmosphere ([Ref. 1]. These localized structures in the 
atmosphere fluctuate rapidly, changing in size, density, location, velocity and 
interaction with neighboring structures. There is also considerable variability 
with diurnal temperature changes, annual weather variations and regional 
climatological! differences. Understanding the mechanisms of this turbulence 
and measuring it on a real-time local level are intrinsic to high power earth 
based lasers. 

These small scale turbulence structures in the atmosphere play an even 
greater role with respect to acoustic absorption and reflection [Ref. 2]. It is well 
Known that acoustic pulses, on the order of 1 to 6 kHz frequency, are ideal in 
defining turbulence structures in the first few kilometers of the atmosphere, with 
lower frequencies having greater ranges. Until recently, the most popular type 
of acoustic sounder was a monopole speaker/receiver mounted in a parabolic 
dish [Refs. 3 and 4]. Several problems with this method include difficulty of 
transportation, high amplitude side lobes at large angles and poor altitude 
resolution due to the low frequencies involved. These deficiencies can be 
overcome by means of planar arrays [Refs. 5 and 6], and this thesis was 


designed to amplify this point. 


Working under Walters, Weingartner and Wroblewski [Refs. 7 and 8] again 
demonstrated the feasibility of this technique by building and testing a 25 
element 5 kHz planar array. The result was an extremely effective low altitude 
echosounder which was less expensive, more portable and had better side 
lobe reduction than previous monople sounders. 

During the course of this previous work several questions arose 
concerning the planar array. If the elements were more densely packed, then 
could the size of the side lobes be significantly reduced, thus increasing the 
accuracy of the beam? Could the array be modeled on a computer to produce 
the best design? Could the enclosure be improved? 

This thesis is an attempt to answer these questions, and involves the 
development of a dense packed 5 kHz hexagonal array and associated sound 
deadened enclosure. Additionally, various computer programs using fast 
Fourier transform algorithms were’ studied as a first approximation in 
diagnosing acoustic side lobe suppression. 

The results indicate that densely packed planar acoustic echosounders 
produce better side lobe suppression than more open arrays and are accurate 


tools in the identification of low level turbulence. 


ll BACKGROUND 


The atmosphere induces random phase _ perturbations on an 
optical beam. Diffraction of these small scale phase variations produces a 
larger beam, restricting the ability to focus energy on a target. Optical phase 
perturbations in the atmosphere stem from temperature fluctuations carried by 
the turbulent velocity field. These temperature fluctuations vary in size, 
cascading from tens of meters in diameter down to centimeters and smaller. 
These structures appear to be most variable and inhomogeneous at the 
interface between stratified layers, particularly near the earth's surface. 

Within the first few hundred meters of the surface, these temperature 
Structures are dominated by the surface heat flux, resulting from differing rates 
of heat transfer between the earth and the atmosphere. During the day the 
Surface heats much more rapidly and non-uniformly than the atmosphere. An 
insulating blanket of air near the surface heats and rises off the ground, more or 
less in sheets. Warmer sections rise more rapidly, quickly disrupting the 
uniformity of these sheets. Cooler air moves in to replace the rising air, further 
disrupting the sheets and causing vertical thermal plumes to form. These 
thermal plumes continue upwards in the atmosphere, slowly mixing and 
increasing the temperature of the lower atmosphere. This cycle is continuous, 
with the actual transfer of energy occurring by small scale molecular diffusion. 


(Ref. 9). 


At night the cycle reverses, with the surface cooling off much more rapidly 
than the air via Planck radiation. In this case, as the air gives up energy to the 
cooler surface, it produces a turbulent boundary layer (via velocity interaction 
with a stationary surface) which remains relatively close to the ground. This 
cooling process Is generally less turbulent than the heating process, due to the 
reduced heat flux at the boundary, although intense turbulence continues to 
exist at the interface of the horizontally stratified layers. 

There is also a neutral event, during the transition between night and day, 
where the surface temperature equals the ambient air temperature. During this 
time the net heat flux between the air and the surface changes sign and passes 
through a minimum. The atmosphere generally becomes optically stable during 
these bref periods. 

The structural variability is clearly evident in optical instruments designed 
to measure the turbulence of the atmosphere [Refs. 10 and 11], but altitude 
correlation with optical instruments alone is difficult. Acoustic echosounders 
alleviate this difficulty. 

Acoustic echosounders are devices which generate and propagate 
energy into the atmosphere and measure the backscatter return from turbulent 


Structures as a function of altitude. The power returned from these turbulent 
Structures, Pp, is proportional to the power transmitted, Py. The equation 


relating these terms has been summarized by Neff (Ref. 12], based upon work 


by Tatarski [Ref 1], and is repeated in equation (1) below. 
(1) Pr= Py ERE? (ct/2) {exp(-28R)} (AG/R2) u[R, f] where, 


- Pp is the electrical power returned from range R, 

- Py is the electrical power transmitted at frequency f, 

- Ep is the efficiency of conversion from acoustic to electric power, 

- Ey is the efficiency of conversion from electric to acoustic power, 

- Cc is the local speed of sound (m/s), 

- tis the acoustic pulse length (s), 

- Ais the antenna area (m2) , 

- Ris the range to the scattering volume (m) 

- exp (-28R) is the round trip power loss due to attenuation where Bis the 
average attenuation to and from the scattering volume, 

- Gis the effective aperture factor for the antenna, 

-u[Rf] is the scattering cross-section per unit volume; that is, the 
fraction of incident power backscattered per unit distance into unit solid 


angle at frequency f. 


The backscattering cross section, yu, is the unKnown quantity of interest in 
Eqn. (1). Tatarski [Ref. 1] uses up in equation (2) (below) to find Cy2, the 
temperature structure parameter. 


(2) u[R,f] = 0.0039 k 1% (Cy2/ T2) where, 
- Kis the acoustic wave number = 2rn/L, 
- Lis the acoustic wavelength, 
- Tis the absolute temperature and 


- C72 is the temperature structure parameter. 
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The temperature structure parameter, C72, is the mean squared 

temperature difference between separate points in space, 

(3) Cre == Malay Myeiyies }/ (Ro- y)22 
and provides the scattering centers which produce backscattered energy from 
the turbulence. <AT>4 empirically increases by AR2%. Therefore, it is normalized 
by AR?* to produce a relatively constant temperature structure parameter. 

Since u is proportional to the backscattered return power, Pr, and to 
the temperature structure parameter, C72, then Pp is proportional to Cr2. 
Solving equation (1) for u and equating equations (1) and (2), the following 
expression for Cy2 is given for an acoustic echosounder of wavelength L; 

(4) Cy2= 1/0.0039 (1/ ERE+) T2 k-'9 (2/ct)(1/AG) Pp/Pr {R2 exp(2BR)}. 

The electric to acoustic efficiency factors for the transducers (speakers) 
were measured in Reference 7 and are inherent to the design of the speakers. 
Ep =~ Ey = 0.5, an excellent value for speaker efficiency, which accounts for 
their use. The aperture area, A, for the 19 element hexagonal array is 0.0866 
mé, 

The effective aperture factor for the array, G, accounts for the non-uniform 
antenna directivity. G is a difficult quantity to calculate or measure. It represents 
a normalized average of the acoustic energy within the transmit/receive 


antenna beam lobe and is clearly summarized by Probert-Jones [Ref. 13]. 


1) 


Based upon Reference 13 and Studies by Hall and Wescott [Ref. 14], a value of 
G = 0.4 is currently used in computations. This quantity is presently undergoing 
extensive research by Dr. Walters. 

Another difficult factor in Equation (4) involves the average atmospheric 
attenuation loss factor, B, which is strongly dependent upon frequency and the 
absolute water vapor in the atmosphere [Refs. 2 and 15]. This dependence is 
illustrated in Figure 1, which is taken directly from Reference 15. The water 
vapor pressure is measured via a digital humidity/temperature indicator, model 


5165-A from Weathermeasure [Appendix A]. 
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Combining these known values into a constant Q, equation (4) becomes: 
(5) Cr2 = (Q) T2k-14 {Pp / (ct)( Py)} R* exp(2BR) where, 
O12 Eeyie00s9 Fy, Ee ET AG). 

Wroblewski [Ref. 8] utilizes this in the software for the echosounder to 
calculate a 15 minute time averaged value for the temperature structure 
parameter as a function of altitude. Uncertainties in the antenna beam shape, 
side lobes, rejection, electronic gains, transducer efficiencies, electronic 
antenna gains and the atmospheric attenuation produce uncertainties in the 
absolute calibration of an echosounder. 

Independent verification of the Cy2 values measured and calculated by 


the echosounder are necessary to confirm the accuracy of this technique. 


Without this confirmation it is still possible to guage relative variations in the low 





level turbulence, simply with the echosounder profiles. These profiles are 


presented and discussed in the following chapters. 


13 


ll. EQUIPMENT 


Much of the equipment used in this project was borrowed _ from 
the 25 element square array designed by Weingartner and 
Wroblewski. The equipment basically involves a Hewlett Packard 
(HP) 217 computer for data generation and analysis, an HP 3314A 
function generator (controlled by the software) which produces the 
initial pulse, an amplifier, a switch/pre-amp which passes the 
initial pulse then amplifies and transfers the return signal to the 
computer, the array and shroud which act as the physical 
transmitter/receiver and a filter used to eliminate stray noise. 
These units will be broken down and individually discussed in this 
section. The software developed for the echosounder by Wroblewski 
was not significantly altered and is not included in this document. 


[Refs. 7 and 8] 


A. ELEMENTS 
The elements of the array consist of 19 piezoelectric 
Motorola brand speakers, model # KSN 1005A, with attached horns 


of model # KSN 1032A. Ceramic piezoelectric drivers were chosen 
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because of their small size, excellent (see Figure 2 ) frequency 
response at 5 kHz, and a shorter 'ring' time with respect to 
electrodynamic speakers. Any Speaker acts aS a damped 
harmonic oscillator with a ‘ring’ time (recovery time) inversely 
proportional to the damping constant. The damping Is much stronger 
in Piezoelectric speakers than In conventional electrodynamic 


elements. 


Typical Frequency Response 






INPUT 28V 
MICROPHONE DISTANCE 1/2 Meter 
{1 Watt into 8 Onms) 





OUTPUT (dB) 


FREQUENCY (Hz) 
200 1000 2000 5000 10000 20000 40000 


Figure 2. Speaker Frequency Response 
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An initial set of 40 speakers was tested to find 19 which 
were similarly matched in frequency response. This testing 
was done by comparing the output of a chosen reference 
speaker with the  acoustic-to-electrical conversion of each 
test speaker. By clamping the reference speaker face-to-face 
with the test speaker, a known voltage to the reference should 
produce a range of voltages from the test specimens. These voltages 
were compared using Lissajous plots on a Nicolet Oscilloscope. It 
was more important to obtain a set of identical speakers rather than 
the best speakers, in order to achieve a more uniform final array. 

The horns on each driver, originally circular with a square 
flange face were trimmed to a circular diameter of 7.62cm such that 
they fitted close packed. Ideally a diameter of d= L/2 
(L=wavelength) would provide optimal phase spacing of elements, 
but this required either smaller speakers (d = 3.4 cm ) which 
were unavailable, or a lower frequency (=2230 Hz), which was 
below the useable range of available speakers. 

Another source of speaker was identified as a Philips brand 


piezoelectric 2002 PT /[Ref. 16] with better frequency response at 
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3kHz than the Motorola speakers in use. Unfortunately these Philips 


Speakers are no longer manufactured. 


B. ARRAY 
The speaker elements were configured in a close-packed 

planar array on an aluminum-covered balsa wood board to absorb 
vibrations (Fig. 3). The elements were wired in parallel, reducing 
the array impedance to 300 2/ 19 = 15.8 Q. The mounting box was 
locally constructed of sheet aluminum and enclosed a layer of 
semi-porous foam loaded with a thin lead sheet to reduce stray 
noise signals. Input/Output was via a coaxial connector. Bracketts 
were attached to allow 360° rotational testing from 2 axes in the 


anechoic chamber. 





Figure 3. Array Dimensions 
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C. ENCLOSURE 

The enclosure was designed to be easily reproducible, 
inexpensive, portable, and easily assembled [the earlier enclosure 
took some time to set up and break down]. Optimally, a long 
cylindrical enclosure was desired to reduce side lobe scatter and 
return without obstructung the main lobe. The enclosure chosen to 
optimize the above requirements was a Rubbermaid 55 gallon plastic 
trash container, with a height of 84 cm, having a slight upward 
taper; the barrel has a lower inside diameter of 53 cm, opening to a 
diameter of 60 cm. This provided a finished open half angle of about 
20°, 8° more than optimally desired according to tests in the 
anechoic chamber with the unencumbered array [see RESULTS]. A 
second possibility was a highway barrier barrel {Appendix A], but 
these were much more expensive and involved similar dimension 
problems. 

The inside surface of the enclosure was covered with lead lined 
semi porous foam, with 7.62 cm convoluted (egg crate shaped) foam 


on one side, and 2 cm flat foam on the backing side. The entire foam 
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liner was attached with standard contact cement. The lead lining 
amounted to a density of about 1 Ib/saft which is sufficient to 
reduce one way noise penetration by 40dB. The convoluted foam 
helped break up internal and off axis scatter (spectral reflection) 
and aided in side lobe reduction. The foam used turned out to be the 
last of a salvage batch of foam, and a second source of foam was 
desired. An ideal substitute for lead lined foam (which was 
difficult to obtain) was found to be a barium impregnated vinyl foam 
(a thin rubbery vinyl layer loaded with barium and sandwiched 
between semi-porous foam) with a similar density of 1 Ib/saft and 
virtually identical sound absorption capability (40dB). Two 


sources of this foam were identified [Appendix A]. 


D. ELECTRONIC HARDWARE 

The data processing system consisted of an HP Model #217 
computer with a 20 megabyte hard drive and 2 megabyte memory, an 
INFOTEK BC 203 Basic compiler, analog to digital (A to D) converter, 
INFOTEK FP 210 floating point accelerator (to enhance the speed of 


execution), HP 9133 floppy disk drive, monitor and HP _ ink-jet 


printer (See Figure 4). The HP 217 was the workhorse of the system, 
triggering the function generator, receiving and processing’ the 
return pulse, and outputting the results in graphical mode. Without 
this unit, timely acquisition and reduction of data would be 
impossible. 

The HP 3314A function generator, adjusted and triggered 
automatically by the HP 21/7 computer, was used to generate a 
Specific transmitted signal. This signal was generally 100 cycles of 
a sinusodial waveform of 1.5 volts (rms) amplitude. 

A QSC model 1700 audio amplifier was used to increase the 
voltage to the array from 1.5 volts (out of HP 3314A) to 30 volts 
(rms). This model was chosen because of its ability to drive a 
highly reactive load (the array). A lighter, more compact amplifier 
could easily replace this item of equipment, and will be considered 
during future upgrades. 

A pre-amplifier (pre-amp) (Fig. 5) was used to isolate the 
transmitted pulse from the received pulse and to amplify the return 


signal by a gain of approximately 11,000. The pre-amplifer was 
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Fig. 4. Echosounder System Set-up. 
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Fig. 5S. Schematic of Pre-Amplifier 
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designed and built by Walters and utilizes two low noise solid state 
operational amplifiers, an LT 1007 and LT 1037. All four isolation 
rectifiers were of type IN 4000. 

The important properties to consider for an effective pre-amp 
were complete isolation of the circuit during transmission (to 
prevent overloading components) and immediate signal retrieval 


following the ‘ring’ time of the array. The more rapid the recovery 
time of the array (a function of driver capacitance (C,) and the loop 
resistance (R,)), the quicker the pre-amp can respond to returned 
signals and the smaller the surface blind zone. For this array, the 
capacitance of each driver was Cp = 0.1 u F. Thus the time constant 
for the array was 
RC = (Cp)(#drivers)(R4) = time constant 
= (0.171078 F)(19)(270 Q) = 0.513 msec. 

For a given array, resistor Ry provides the time constant of the 


Sire uit 
After 5 time constants there is virtually no ‘ringing’ left on the 


array, and this sets the minimum height of the blind zone. Thus: 
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5xT=2.57ms and Happ = (340 m/s)(5.32710°3s) = 0.87 meters. 


Therefore the blind zone stands at approximately 1 meter above the 
array, which is not noticeable on the echosounder printout. 

A Rockland Wavetek model #753A filter {115 dB/octave} was 
used to isolate noise and further amplify the return signal prior to 
computer analysis. The input side of the filter was set up as a low 
pass filter at 5 kHz and gain setting of O dB (x1); the output side 
was set up as a high pass filter at 5 KHz with a gain setting of 20 dB 
(x10). This produced a band pass filter centered at 5 kHz with a 
band width of about 100 Hz at the -3 dB level. The over all gain of 
the filter in this mode was measured as 

Output 685 m volt 
ed = «9G = 10 
input 70.88 m volt 
The combination of the pre-amp and the filter produced an overall 
return signal gain of (pre-amp gain)(filter gain) = (11094)(9.67) 
= 107,000. 
The gain of the system is thus: 107,000 * 270 Q/( 270 + 15.8 Q ) 


= 101,000. 
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IV. RESULTS 


A. LOBE PATTERNS 

The array was tested in the anechoic chamber utilizing a 
program designed by Butler [Ref. 17]. This program utilizes the 
set-up shown in Figure 6 and allows the measurement of beam 
patterns for various arrays. A band pass filter was used to increase 
gain and reduce 60 Hz noise from the equipment. Present on all the 
lobe patterns is a spurious lobe at = 135°. This was determined to 
be a reflection from the inside door handle of the anechoic chamber 
and should be ignored. All lobe patterns have been normalized to the 
maximum voltage detected, allowing easy comparison between 
figures. 

The hexagonal array was rotated around two axes to obtain the 
maximum sidelobe field of the array (Figs. 7 and 8). The pattern for 
the square array is shown for only one axis of orientation (Fig. 9). 
Comparison between the two arrays (Figs. 7 and 9) shows a definite 
decrease in side lobe number and intensity for the hexagonal array 


versus the square array. 
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Figure 6. Anechoic Chamber Set-up 
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While the array was being tested in the anechoic chamber a 
demonstration of its detecting ability was desired. Using a heat 
gun, the array was set up in the chamber exactly as it would be in 
the field, except a Nicolet model #3091 digital Oscilloscope was 
used to detect the return signal. 

The heat gun was tested at three different settings; first was 
‘heat and air’, then ‘air only’ and last ‘heat only’. With the heat and 
air both on, a definite fluctuating signal could be detected on the 


oscilloscope screen. With ‘air only’ there was virtually no 
detectable signal and with ‘heat only’ the signal was almost back to 
its level of ‘heat and air. No measurements were taken, but the 
demonstration was repeated several times. This directly illustrated 
that the turbulence structures in the atmosphere are mainly a result 
of temperature fluctuations and not due- to turbulent velocity 
fluctuations. 

The next step involved testing the side lobe suppression of the 
enshrouded array. The mechanics of this operation were 


considerably more involved than simply suspending the array from a 


rotating motor. Being too heavy to suspend, the shroud was taped to 
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a rotating turntable and centered beneath the rotator motor. The 
motor was connected to the shroud via aluminum rods and clamps. 
When the rotating motor turned, so did the shroud, allowing the 
enshrouded beam pattern to be measured and recorded with the same 
software as used previously. 

The initial beam pattern, (Fig. 10) showed that the shroud, while 
significantly reducing side lobes, was not long enough to completely 
eliminate them. Some of the side lobes were getting through. 

A first-cut step in further reducing these side lobes involved 
placing a 15 cm wide sleeve of barium impregnated smooth foam 
around the upper inside circumference of the shroud. The result (Fig. 
11) was not promising, and demonstrated that further refinement 
was required. One thing Figure 11 demonstrated, however, was that 
specular reflection from a smooth surface increases side lobe 
activity and thus convoluted foam shrouds (vice smooth foam 
shrouds) are required. 

The next step involved placing a 15 cm wide annular ring coated 
with sound deadening foam over the exit of the shroud (similar to a 


trash can lid with the center removed). This left an open circle of 
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48.25 cm diameter for the acoustic waves to propagate. This 
adjustment slightly decreased the side lobes (= 3 dB), but not 
significantly (Fig. 12). Further narrowing of the neck via more 
annular rings only succeded in reducing the main lobe, with 


little subsequent change in side lobe activity. 


B. FIELDMEASUREMENTS 

Field measurements were conducted on several occassions 
over a six month period at various locations around California. They 
generally coincided with optical measuremenis being conducted by 
Walters [Refs. 9 and 10] and were used to supplement the optical 
data. Of greatest interest were mountaintops above the marine 
boundary layer (=4000 ft.), where the turbulence should be a 
minimum. 

Two mountaintops in particular were visited, one being Lick 
Observatory on Mt. Hamilton and the other being Mt. Wilson 
Observatory at Mt. Wilson. These were chosen because of the 
availability of telescope domes, which were required to reduce wind 


vibration on the optical components and protect the computers from 
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dust and debris. The echosounder itself is impervious to wind, and 
can be operated wherever electricity is available. 

Lick Observatory was visited on 1 and 2 August 1987 and some 
of the results are presented in the following figures. All times are 
given as universal time. Subtract seven hours to obtain local 
(daylight savings) time at the site. Both days were similar with 
respect to turbulence. Therefore, only profiles from 2 August will 
be presented to illustrate the ability of the echosounder to detect 


and present the rapidly varying nature of the atmosphere accurately. 
The C72 profiles are currently unverified and should be ignored. 


Figure 13 was taken early in the evening and demonstrates the 
presence of convective pluming. This pluming vanishes gradually, 
eventually yielding to the nighttime structure which starts in Figure 
14. By this time a definite layer has set in and surface pluming is 
no longer dominant. Later in the evening a strong turbulence layer 
descendes (Fig. 15) which quickly settles [Fig. 16] and is quite dense. 
This layer remains close to the ground the rest of the evening, with 


light turbulence detectable above it. 
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Fig. 13. Echosounder Profile; 1830 - 1845 Local Time 
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Fig. 16. Echosounder Profile; 2215 - 2230 Local Time 
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Early in the morning (Fig. 17) another strong layer is noted. As 
the sun begins to rise (Fig. 18), commencement of the daily 
convective pluming becomes clearly evident. These figures 
illustrate the dramatic ability of the high frequency echosounder to 
profile low level turbulence in the atmosphere. Interestingly, 
whenever a turbulence structure was detected on the profiler, there 
was also noticeable degradation to both the optical instruments in 
use. 

The next occasion to use the echosounder occurred on 6 and 7 
November 1987 with a field trip to Mt. Wilson Observatory. 
Unfortunatley, no echosounder data could be obtained due to severe 
radio frequency [RF] interference from approximately 10 radio/TV 
transmitters co-located at the observatory. The trip was not a loss, 
however, as it was used to correct a subtle flaw in the design of the 
pre-amp. 

The problem was traced to the ‘floating point ground’ design 
of the pre-amp. The coaxial cables leading to the pre-amp made 
ideal antennas at the 200 megahertz range, transferring enough 


radio frequency voltage into the pre-amp casing to influence the 
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circuit. The rectified voltage effectively overloaded the first stage 
amplifier, overriding the weak return signal from the array and 
destroying the data. This effect went undetected in previous use 
of the pre-amp. The problem was corrected by improving the 
grounding structure of the device. 

Note that all these profiles were taken with the hexagonal 
array mounted in the square shroud, as the hexagonal shroud was not 
yet completed. To ensure that the quality of data was not 
compromised, the square array was also run at the site and the 


profiles were identical. 
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V. ANTENNA BEAM PATTERN MODEL 


This thesis also modeled the planar acoustic array as an 
optical aperture, obtaining the Fourier transform of the far field 
‘diffraction’ pattern via computer calculations. The technique rests 
on the fact that 

"". , the Fraunhofer diffraction pattern is the Fourier transform 

of the field distribution across the aperture... .". [Ref. 18, p. 412] 

The quote applies to electromagnetic radiation, but Dr. walters 
correctly construed the electromagnetic amplitude as an acoustic 


pressure, obtaining the following equation: 


(6) — P(k,,k,) = 1/iR JJ P(x,y) exp {i (k,x + ky)} dxdy ; 


aperture 

- P(ky, ky) is the pressure distribution in the image plane, 

-k, and ky are the spatial frequencies along the x and y axes of 
the aperture respectively; kj= k*J/R, where J= x or y; k=2n/A , 


- P(x,y) is the pressure amplitude at the aperture, 


- exp {i (kK,x kyy)} is the spatial variation of the aperture and 


- R is the distance to the image plane; | accounts for temporal 


phase distributions across the aperture. 
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The model can be thought of as a plane coherent acoustic wave 
passing through the aperture, which produces a diffraction pattern 
(the Fourier transform) in the far field (the image plane) by means 
of the Huygens-Fresnel principle. From this pattern it was possible 
to reasonably predict the size of the main lobe and the approximate 
angles and intensities of the side lobes. This modeling was done for 
both the square array and the more densely packed hexagonal array. 


The results are presented in this chapter. 


A. HARDWARE 

The computer used in this phase was a Hewlett Packard (HP) 
217 with an INFOTEK Basic compiler, identical to the computer used 
in the echosounder arrangement. The memory requirements for this 
program, however, were tremendous, requiring over 4 megabytes of 
memory when using the 512 x 512 grids. This lead to insertion of 2 
additional memory boards into the computer, (one from INFOTEK and 


one from EVENTIDE), bringing the hard memory to over 5 megabytes. 
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B. SOFTWARE 

The program [Appendix B] utilizes a Fast Fourier Transform 
(FFT) routine supplied by INFOTEK. The idea is to enter the aperture 
(array) aS a small series of points on a huge grid. Typically the grid 
size would be 256 x 256 points square, and occassionally 512 x 512. 
The array takes up only a small center portion of this grid, 
approximately 40 x 40 grid points. The larger the grid field in 
relation to the input aperture, the better the resolution, provided the 
input aperture is not so small as to obscure the finer details. 

The aperture area was entered digitally as a series of ‘ones’, 
Signifying a uniform coherent plane wave at the aperture, with the 
remainder of the field entered as zeroes. The FFI routine then 
computed the far field Fourier transform of the array, and plotted 
the output intensity in relation to a gray scale. The gray scale was 


adjusted to encompass 3 decades of intensity. 


C. RESULTS AND ANALYSIS 
A sample of the input apertures, both square and 


hexagonal, are shown in Figures 19 and 20. Generally, the radius 
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of each element is only 4 grid units, with a background field of 256 
x 256 grid units. These figures have been enlarged to show detail. 
The output Fast Fourier Transform for each is shown in Figures 21 
and 22 respectively. The center portion of the FFI for each array is 
the lowest spatial frequency for the array, and corresponds to 
the main lobe. Rings further out correspond to the side lobes at 
larger and larger angles. Rings far from the main lobe are not 
physical except at very high frequencies. Interestingly, as the 
frequency is increased, more and more side lobes are forced into the 
forward hemisphere. 

The program calculates the beam structure of a single piston 
driver on an infinite plane. The actual source is not a piston driver 
(but the approximation is close). Huygens wavelettes propagating 
backwards from the end of the horns reflect from the support plane, 
increasing the energy of the actual lobe field. Also, the physical 
side lobes can extend no further than 90° from the axis of the array 

(due to surface absorption and reflection) but this limitation is 


ignored by the program. In addition, an obliquity factor was not 
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included. These limitations affect the graphical output, but do not 
significantly affect the following calculations. 
The angles of the output beam pattern rings depend on the 


frequency of interest (5 kHz), the dimensions of the actual array and 


the input/output dimension of the aperture. To calculate 


information from these printouts it is necessary to use optical 


aperture equations from Hecht and Zajac [Ref. 18]. The half angle 


formula for a circular aperture is: 


(7) sin@ = Q/R = 1.22 L/2a where, 
- q is the distance from axis to side lobe of interest, 
- ais the radius of aperture, 
- Lis the wavelength, 
- R is the distance to the far field, 


- © is the half angle or cone angle from the axis to q. 
The hexagonal array is not quite circular, but in the limit it 


approaches a circle. In this case the constant 1.22 is simply 


replaced by F, a constant to be determined by experiment. 
By using Brigham's text [Ref. 19] it is possible to convert 


equation (7) into the form: 


(8) O© =F ARCSIN (Ln/NAT) where, 
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- @ is the angular measure of the side lobe from the main lobe, 
- Lis the wavelength of interest = 6.86 cm, 
- NAT is the total width of the sampled grid = (ND/X), 
- N is the width of the input window (grid units;128,256 or 512), 
- Dis the actual dimension of the physical array (cm), 
- X is the dimension of the input array in grid units, 
- n is the output number of grid units to the side lobe of 
interest=(rz/w), 
- r is the distance to the output side lobe of interest (cm), 
{NOTE: X and D must correspond; r and n must correspond}, 
- z is the width of the output window in grid units (< N/2), 


- w is the width of the input and output windows in centimeters. 
As an example, consider Figure 23. For this example, D = 
38.735 cm (measured between the outer edges of the farthest 2 
speakers), w = 10.3 cm, N = 256 grid points and z = 128 grid points. 
To calculate the angular extent of the main lobe in real space, first 
measure the distance (in cm) from the corner of the output quadrant 
to the edge of the main lobe, along the horizontal axis (Fig. 23). 


Here, r = 0./ cm and X = 40 grid points. Therefore, equation (8) 


becomes: 


(6.86) (0.7)(128/10.3) cm 
(9) © = F x ARCSIN ----------------------- =| oad am oS ae 
(256)(38.735)/40 cm 
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Figure 23. Output Quadrant of Hexagonal FFT 


The half-angle of the actual array is only 12° as measured on 
Figure 8. This means that F = 0.86 or 1/F = 1.16 which its very close 
1On tucc. Using the value F = 1.16, then the 90° point can be 
calculated for the output field. In this case, 

(2565)(38.735) 


(10) MUO == =e" See Gi 
(6.86) (128/10.3)(40) 


3)! 


which means everything outside the corner radius of 2.81 cm can be 
ignoredas it is non physical. 

The output FFT'’s were correlated to the actual lobe patterns of 
Figures 7 and 8, the apex and side mounted patterns for the 
hexagonal array, and the results are plotted in Figures 24 and 25 
and tabulated in Table 1. The amplitudes were normalized by 
assigning a value of unity to the main lobe; the angles are measured 
counter-clockwise from the standard X axis; the amplitude scale is 
in decibel units. 

The measured amplitude was calculated using INVLOG (dB/20), 
where ‘dB’ is the normalized decibel level estimated from Figures 7 
and 8. The computed amplitudes relied on the formula 20 LOG 
"Q/Max" [see Appendix B, lines 1390 and 2380 for definitions of Q 
and Max]. 

The angles were computed using a simpler form of equation (8); 

(11) O = Arcsin { (K-1) L/NAt }, where K is an index from 
1 to N/2, and essentially allows the computation of the angle at 
each grid point chosen. The measured angles were taken directly 


from Figures 7 and 8. 
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Fig. 25. Polar Plot of Side Mounted Hexagonal Array 
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TABLE 1. COMPARISONS BETWEEN ACTUAL LOBE PATTERNS AND 
COMPUTED LOBE PATTERNS. 


APEX MOUNT (Fig. 7 vs. 24) 





Measured Computed 
20a amplitude © __amplitude 
1St min ting 025 12s 002 
1St max Vee 126 1S 128 
2nd min 21° 056 23° .007 
2nd max 26° 095 28° 108 
3° min 33° >.01 38° .006 
3°d max 42° 056 41° O27 
4th min 51° 018 46° 007 
4th max 64° 056 62° 148 


SIDE MOUNT ( Fig. 8 vs. 25) 


Measured Computed 
Q_ amplitude a amplitude 

1St min Vz .050 14° 048 
1St max 18° 158 19° 151 
2nd min Dig >.01 28° 01 

2nd max 35° 039 36° 030 
3d min 52° >.01 46° 024 
3°d max 61° 025 56° 029 


Ole 


The side mounted section corresponded to measurements along 
the X axis, and the apex to measurements along the Y axis. The error 
of the measurements was grossly estimated as approximately 10% 
for the minimums, as they tended to be obscured, and 3% or less for 
the maximums. Some values exceed this amount considerably and 
require further investigation. 

Generally, however, the tabulated values agree strongly and 
Suggest that this technique has considerable credibility at 


predicting side lobe structures of untested arrays. 
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VI. CONCLUSIONS 
This thesis demonstrated that increasing acoustic planar array 
Speaker density from 25 speakers/1452 cm? (1.72 x 10°° 
speakers/cm2) to 19 speakers/943 cm? (2.02 x 10°° speakers/cm?) 
effectively Suppressed the lobe spreading from a 70° cone angle to 
less than a 45° cone angle, with a reduction in side lobe intensity of 
between 4 and 10 GB. 

Graphical profiles demonstrated the great variability and 
volatility of turbulence structures in the atmosphere and underline 
the need for better tools to measure and understand this structure. 
These profiles identify gravity waves, thermal plumes, boundary 
layers and other temperature induced structures. Particulary 
important, this device allowed tentative quantification of altitude, 
thickness and relative intensity of mountain boundary layers over a 


potential DOD laser site. 
Independent verification of the Cy* values is important and is 


currently under study. It will further enhance the utility of this 


device. Due to the strong exponential dependence of acoustic 
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attenuation (8) with increasing frequency, future improvements 
Should also include lower frequency drivers (2-3 kHz) of similarly 
high efficiency, to improve the range of the device. 

It was shown that acoustic array beam patterns could be 
successfully modeled on a computer by means of a Fast Fourier 
Transform algorithm and tested for optimal design. Future 
improvements to this technique include a method to stochastically 
simulate turbulence and model the beam structure as it passes 
through this artificial turbulence. Calculation of main lobe energy 


versus total energy for various arrays should also be studied. 
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APPENDIX A 
List of Material Sources 
Disclaimer: These sources are listed only in case exact 
duplication of this thesis is required at some future date. Other 
equally competitive sources for each item do exist and their 
exclusion is not a reflection of product quality in any way. 


Barium impregnated foam: 


1. Barrier Corp. 2. Latchford Corp. 
Tigard, Oregon (714) 734-3660 
(503) 639-4192 POC: Tom Pelligreeni 


POC: Mark Dove 


Cylindrical Enclosure: 


1. Rubbermaid, Inc. 2. Traffic Products Co. 
1147 Akron Rd Citrus Heights, Ca 
Wooster, Ohio 44691 (916) 723-2958 
(216) 264-6464 POC: Michael Fish 

Speakers: 


Motorola Inc. 

4800 Alameda Blvd N.E. 
Albuquerque, N. M. 87113 
(505) 822-8801 


Temperature/humidity !ndicator 
Weathermeasure 
P.O. Box 41039 
Sacramento, CA 95841 
(916) 481-7565 
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APPENDIX B 


Fast Fourier Transform Program 


12 | Program plots either the input points, the FFT or the 
Polar Plots of a 19 element hexagon array or a 25 

Aa | element square array. 

30 | =6This program 1s interactive. 

4 Pi Wec., 16. EE/LRM modified in OCT 87 

5 momen Cl. oie) imc | 351 Z) 

60 PovecerR fan 


73 Cee d ese S12 ).imag Siz 5l2) 
BQ REAL Amp( 100) ,Ang(10Q) 


1) | 

1Q2 Boece Eee GRID SIZE N; N= 5l2,256 OR 128" 
11@ INPUT N 

12Q PRINT 

130 | 


140 NE=N/2 

15@ Mal OGrwN)/_OGG2Z) 

16@ | 

17@ foie ENIER * 1’ 10 SEE THE 19 ELEMENT HEARGON. ” 
182 pee fee NtER “2° TQ SEE THE 2h ELEMENT SQUARE." 
190 INPUT xX 


2QQ Peng 
210 PRINT 
a. re eaten TO eSEE A PLO! OF THE INPUT POINTS" 


2 58 teens 2 lO SEE A PLOT OF THE BRT” 

240 Poe wee  Saeio SEE A POLAR PLOT Sr TENSITY" 
Bod INPUT A 

26@ PRINT 

2(Q PRINT 

282 Pode iene GRE SIZE OF EACH BLEMENT, R.~ 
a0 geeiee USUALLY R=4 OR 8.” 

32Q INPUT R 

31 PRINT 

oe PRINT 

5.50 apie ENTER? 7 TO SEE THE FIRST QUADRANT ONLY” 
548 are ecloe ENTER *@°.” 

359 Nee 16 


362 Mme-esie SsiZe OF eIHE OUTPUT WINDOW 
378 B=N2/2 

588 | 

398 [F X=1 THEN 

402 | 


412 poe R* COSC. 5255988) 
429 FOR I=N2-6*R TO N2+G*R 


430 FOR J=N2-G*R TO N2t+6*#R 
44Q IF ((1+4#R-N2 )°2+(5-N2)°2)<R*R THEN 
4590 Real( I ,J)=1.000 
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460 
47@ 
480 
49 
500 
51@ 
520 
530 
540 
S50 
56 
57) 
580 
S30 
60Q@ 
610 
B20 
650 
640 
650 
660 
670 
680 
69@ 
7100 
110 
720 
750 
740 
750 
760 
7710 
782 
792 
88@ 
819 
820 
830 
840 
850 
860 
870 
830 
89@ 
98@ 
ao 
920 
950 
94 
95@ 
960 


FOR K=)) eee 
Real(I+2*kK*R ,J)=Real( I,J) 
Real( 1+2*h*kR-Re Ith =Realt i ai 
Real(I+2+*#K*R-R ,J-F )=Real( lI J) 
NEXT K 
FOR GK= 20 
Real( 1l+2*kK#R -J=246 D=Reali 1 J) 
Real(I+2*kK#R J+2*F )=Real( I,J?) 
NEXT K 
END IF 
. Ge 
NEXT od 
| 


END IF 
| 
IF X= 2S THEN 
| 
K=N2+54#R 
a2 Sek 
EO Ra— ie 0 a 
FOR J=| TOUK 
IF I1<K AND I2L AND J<RK ANDO JOL THEN 
Cy=N2-4+#R 
POR P= let.0- S 
Cx=N2-4+#R 
POURS O—=| 10°55 
IF I<Cx+R AND I>C.~-R AND J>Cy-R AND J.Cy+R THEN 
[Rea eR OU Cl-0. 2 einen 
Real(J,1)=1.000 
ee 30 
END IF 
Cxa=Cat2e*R 
NEXT O 
Cy=Cy+2+#R 
. NEXT = 
ENG 
Eka) 
NEAT 1 
END IF 
PRINTER MiSs 3701 
| 
| This section plots the input points only 
| 


IF A=1i THEN GOTO 2242 
| 
|. THIS@SECILON 2S ae sr aol POUT east O) min 
| 
TI=TIMEDATE 
FOR [=1 TON 
FOR J=i TO N 
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970 

980 

990 

1002 
1Q10 
1@2@ 
10352 
1042 
10S 
106V 
1070 
1@eo 
1092 
110u 
1110 
1120 
| 3 
1140 
1359 
116@ 
1170 
1180 
1190 
1200 
1210 
lee 

1250 
1c4@ 
1258 
1260 
12 7@ 
2380 
he 38 
1500 
| Sa 
Seo 
5!) 
1340 
Lj3S@ 
1362 
137@ 
1580 
1590 
1400 
1412 
1420 
14352 
1442 
1450 
1462 
1472 


Re(J)=Real(I,J) 
Im( J )=Imag(I,J) 
NEXT J 


Gaeeer rt (Mm) 


POR J=1e'O -N 
Real(I ,J)=Re(J) 
baat lJ =Imi J) 

NEXT J 

NEXT I 
| 
ele =i) Tele 

FOR I[=1 TON 
Re( I )=Real( I,J) 
bit =lTmagt 1, J) 

eel 


Coe ee eh 


POR f=(lF a0) N 
Real( 1 ,J)=Re(I) 
Fmagtl J d=iIm( 1) 
Nesp el 
NEXT J 
ie MEOATE 
Pee INICIO 12-71 )7b0)/100, “MINUTES” 


max=(1+INT(D*N/(68.6#R))) 


| 

Meee - THEN GOTO 1960 

| 

(This section graphs the polar plot of the FFT. 
| 

BI=N2/2 

B2=] 

moeeimeractuad dimension of the @arraytcm) 

we. (55 

P 


The 'K’ loop does the side and apex mounts 


Ri=Real(!,1) 
la=Imag(l ,1) 
Maa=SQR(R1+*Rit+latla) 
| 

0 5 ee) Ones 

| 


IF A=! AND K=1 THEN 


go Uh HEXAGONAL» ARRAY, SIDE MOUNTED" 
ENOSIE 
Hee A= AND K=2c THEN 

pi ae Be neONoiesneRet,, MFA MOUNTED” 
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1430 "END Sr 
1480 IF X=2 AND = eer, 


1580 Seo SQUPRE ARRAY = STOE TOUNTER- 
(S10" “ENG e 

1Ss20. ! 

1530 ! This section computes the ape« square array 
aco 

1550 TF K=2 ANOPK =2e aed 

1560 PRINT SQUARE ARRAY. APEX MOUNTED 
1578 ! The ‘L'toop does only the apex mounted square 

1580 FOR =] 1GNe 

[S30 Rl=Real (ere) 

16G8 6G=Imag(L,L) 

1612 Q=SQR(R1L*R1+6*G) 

1622 IF L>Pmaa THEN GOTO 1660 

1654 Amp(L)=Q/Maa 

1640 Ane (i =ASNK( (i= 19*6855"ho eNO eo PI 

1650 =| PRINT _ ,Q/Ma«,ASN( ((L-1 )*60.6*R 0/7(D2N))*160737 14 1Soa 
1660 NEXT L 


1670. CALBePolarag ly oanat* Ang «ema, 
1682 TE X=. THEN SOTO. 2540 

1630 END IF 

1700 | 

I7res 

[728 FOR P= yi Geet 

1730 POR? J= rows 

1740 IF Bl=|> FHENSP 

1752 IF @2=1. SHEE =) 


1762 Rl=Real(I,J) 

1770 G=Imagi I,J) 

1780 Q=SQR(R1#R1+6#6) 

bv a0 [> PHEMenw THEN GeTe 1640 

1822 Amp(P )=Q/Max 

1812 Angt( FP Y=ASNC CC P=1)*68-6*Ro7 (ie) 207 

is2ge | PRINT P ,Q/Max ,ASN(((P=1 0*68 6*R07(0*N) 8180/7 sn 1415 sae 
1830 NEX Td 

184@ NEXT I 

Reis), ae! 

1860, CALL Polaralt came (* seme or mae 

1870 | 

188@ B82=B1 

189@ Bt=1 

1900 NEXT K 

131@ 

192@ IF A=3 THEN GOTO 2590 

1330} 

1948 ! This section centers the FFT onto the output grid. 
l350% ot 


1360 IF C=7 THEN ! USE ONLY TO GE? ONE QUADRANT Obes HE vcr 1p 
197 GCLEAR 
1980 GINIT 
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J 
J 

} 
S 


Pom Bl POOR Fo fo BD Po oN a 
fl Gh Gl Ol Gre Cr mori hy taf fo EF 
CAE Gl es 1 Oe) a ee 
Soaqoe e_eragandoaagn.seeosa sS 


GRAPHICS ON 
VIEWPORT 30,130,0,102 
WINDOW 


B1l=Ne-1 


oe 
GOTO 


END IF 


2340 


la 


| NZ 


IRS section centers the output. 


pemer=1 10 N2-t 


FOR J=1 


TO Né-1 


Rl=Real( I,J) 
la-lmagi ].J> 


Reo Ne tie N2+J—-|J=R1 
Real(N2-I .N2+J-1)=R1 
Real( N2=-1 ,N2-J)=R1 
meceChear f=  N2—J )=kR | 
Imag(Ne-I ,N2-J)=Ia 
magi Ner l=) Ne-d )=la 
imagenctl— | N2tJ—lj=la 
Enaethe-loN2t+J=|)=la 


NEXT 


lext 1 


GCLEAR 
SIL 
GRAPHICS ON 


{ 


ViEewPOR| 30,152 ,2,10@ 
WINDOW NoO-B N2Z2+8 ,N2-8 ,N2+B 


Pape Tene) TO Note? 
eae J=N2-B 1 


J 


TOSne+e2 


Rl=Real (I,J) 
Ta=Imag( I,J) 


Q=SQR(R1I#*Ritla*la) 
THEN 
Pe( 20s) yy Vs 


iF 


A= | 


END IF 


te 


A=¢ 


END IF 


eee sOSfHEN 2=0 
THEN Z=1 
GeemeONIENSITY 2,2 ,2 
Move 1 , J 

RECTANGLE 


ie 


oo 


THEN 
Z=|+L6T(Q/Max« )/D 


lrdlee eeLee 


| 


Use ONE ORs e B@merHe 


Pe OTeeOlTs 


CB SEORSDSEADE INTENSITY SPREnD 
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2500 NEXT J 

2510 NEXT I 

Zaeoe 4! 

2530 PRINT “The size of the Outpur Windows soe cee 
2540 PRINT “The size of the input window 15° .N."X” N 
2SSQ PRINT “The radius of each element is ",R 

2560 =! 

2570 ! Pause statements allow insertion of parameters w/o 
2580 ! redoing thesFrl calem ation. 

25390) PRINTER Toe) 

2600 PAUSE 


Sei. FPiause 
eoewe Ep 
c650 =! 


2640 #! This section 15 the FFT subroutine (icomeinr ole.) 
efou 

Seo0 SUBPE tthe Gerais) 

2570 COM Rel i:St2 ye inOi: sie. 


-68@ | FET SUBROUTING FROM INEOTEK we seine en coe 
2690 | USE COMMON FOR FASTER EXECUTION s. UsMscecm ces. so 
2700 INTEGER 1° JK Net pee te 
aan’) N=2°M 

2712Q J=] 

2750 POR 1=1 TQ Ne] 

2740 beeen 

Lee T=Re( J) 

-76@ Re(J)=Re( 1) 

ctf Re(I )=T 

a1eU T=Im(J) 

2790 Im( J )=Im¢ 1) 

2800 Im( I )=T 

ae 09) ENGae 

2B2Q K=N DIV 2 

2830 Wine 

284@ esas 

£2850 K=K DIV 2 

2860 END WHILE 

28790 J=J+K 

28B8Q NE Ae leel 

~839@ Le=| 

2908 Olas || POL 

2910 Lel=Le 

2920 Le=LetLe 

Fags) OG) Ure=] 

2940 Uim=2 

2950 Ang=PI/Lel 

29608 Wre=COS(Ang) 

2970 Wim=SIN(CAng) 

298@ FOR J=1 TO bet 

23990 POR lad RON Glee ie 
500 Ip=I+Lel 
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J] 
ud 


Cl Chay ta rat ot) Poe 


1 GO) GF Oo! OG) GF Ol GW oF Gl GI ol Gl ca cad Gl OI) 


JO) GO) Gl GO! Ga UF 


teeehke  ip)*Ure-Im( Ip )*Uim 
Tim=Re([p)*Uim+Im( Ip )*Ure 
Re( Ip )=Re(I)-Tre 
Peel )=lm( 1 iim 
Re(I )=Re(I)+Tre 
ale) = Tia fe iam 
NEXT I 
T=Ure*Wre-UimeWim 
Uim=Ure*#Wimtlime*Wre 
Ure=T 
NEAT J 
Meme 
Noscale=1./N 
MAT Re= (Npscale)*Re 
NAT Im= (Npscale)+Im 
SUEEND 


Mase secvtom 25 the polar plot subroutine written 
pyoreul Uavison. 


_— —_ -_— 


suber olaratt (Amp *) Ang(*),Pmax) 

REM FF FtHtt tee te tet te te tetetteettegtt+etetrtets+ ete 

REM PibometUee NWNArS IHE CALCULATED, NORMALE ZED 

REM PSstseany Dale) ON BOG NORMAL POLAR PLOTS. 

REM FtHttHt Het ttt tt ete tate ttetetretetettteteee tet 

eet ON | ieee NG TRE PLOVER t++4++++4+4+4++++44++ 

GINIT 

Pmeom@ter lo CRT. INTERNAL’ 

GRAPRECS ON 

St Penk 

| ++4++ SCALING THE GRAPH +++t+t++t+ttte+t¢ttet+tetett 
| 

VEEWPORT 25,100 ,20e9S 

SHOW -40 40 ,-4@ ,40 


mile, trre 4 
Ring=!10 
picomhing Ring,” ,0,5,5,2 
Siete. 3 
DEG 
BORG 7 
ColZe 3 
FOR R=@ TO 4 
Sheer Yre | 
MOVE @,R*Ring 
SpSeerk thi ng-40 
Meer TYPE 3 
FOR Angle=@ TO 360 
Pee ieee hime *COS(Anole) R*Ring*SINi Angle) 
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3520 
5530 
3540 
S550 
3560 
5570 
5580 
Saa0 
5600 
3610 
5629 
5639 
5640 
3650 
5660 
3672 
568 
5692 
3700 
3712 
5720 
5750 
3740 
3750 
3762 
3772 
3780 
S740 
5800 
5812 
552Q 
5830 
5840 
565@ 
5860 
387 
53880 
5830 
3938@ 
591 
3920 
5330 
39340 
J330 


NEXT Angle 
PENNE 
NEXT R 
FOR Angle=30 TO 15@ STEP 39 
IF Angle<x>930 THEN 
PLOT 40*COS( Angle) ,40*SIN( Angle) ,-2 
DRAW 40+*C0S(180+Angle) ,40+*SIN( 18@tAngle) 
Ee Nii 
NEXT Angle 
PENGIE 
| 
Eerie 
DEG 
FOUR P=Pmax 10°) sheer —1 
R=40+2Q0*LGT(Amp(P )) 
IF R>=0 THEN 
PLOT R*COS( Ang (PO ResINCAne ce) 
EIN ae ites 
NEXT FP 
FOR P=). FCs erie. 
R=40+20*LGT(Amp(P)) 
[PRS Oe tne 
PLOT R*€0S( Ang. P o> -ReSiik Aner? 
ENO tei 
NEXT P 
| 
CEIP Gee 
l 
LORG 2 
MOVE 40,0 
EAB Els xa 
LORG 8 
MOVE Ago 
PAREN 2 aa 
EGRG6 
MOVE @,-4@ 
LABELER = %- 
| 


CETeeon 

| 

DUMP SOR Ar nh UGoe wo 
GRAPHIC s ORE 
PRINT CHRS( 12) 
SUBEND 
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